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We show that polar molecules driven by microwave fields give naturally rise to strong three-body 
interactions, while the two-particle interaction can be independently controlled and even switched 
off. The derivation of these effective interaction potentials is based on a microscopic understanding 
of the underlying molecular physics, and follows from a well controlled and systematic expansion into 
many-body interaction terms. For molecules trapped in an optical lattice, we show that these inter- 
action potentials give rise to Hubbard models with strong nearest-neighbor two-body and three-body 
interaction. As an illustration, we study the one-dimensional Bose-Hubbard model with dominant 
three-body interaction and derive its phase diagram. 



Fundamental interactions between particles, such as 
the Coulomb law, involves pairs of particles, and our un- 
derstanding of the plethora of phenomena in condensed 
matter physics rests on models involving effective two- 
body interactions. On the other hand, exotic quantum 
phases, such as topological phases or spin liquids, are of- 
ten identified as ground states of Hamiltonians with three 
or more body terms. While the study of these phases and 
properties of their excitations is presently one of the most 
exciting developments in theoretical condensed matter 
physics, it is difficult to identify real physical systems 
exhibiting such properties - a noticeable exception being 
the Fractional Quantum Hall effect. Here we show that 
polar molecules in optical lattices driven by microwave 
fields give naturally rise to Hubbard models with strong 
nearest-neighbor three-body interactions, while the two- 
body terms can be tuned (even switched off) with exter- 
nal fields. 

The many-body Hamiltonians underlying condensed 
matter physics are derived within an effective low en- 
ergy theory, obtained by integrating out the high energy 
excitations, fn general, this gives rise to interaction terms 



K« ({r j) - J2 V (r i -r j )+ £ W (r,,^, r fe )- 

i<j i<j<k 



(1) 



where V(r) describes the two-particle interaction de- 
pending only on the separation between the particles. 
The second term W(ri, Tj, rfc) is the three-body inter- 
action, which depends on the distance and orientation 
of three particles, and vanishes if one particle is far 
apart from the other two. The ellipsis denotes possible 
higher many-body term terms. While for Helium atoms 
in the context of superfluidity the two-particle interac- 
tion dominates and determines the ground state prop- 
erties with the three-body interactions providing small 
corrections,^ model Hamiltonians with strong three-body 
interactions have attracted a lot of interest in the search 
for microscopic Hamiltonians exhibiting exotic ground 
state properties. Well known examples are the frac- 
tional quantum Hall states described by the Pfaffian wave 
functions which appear as groun d sta tes of a Hamilto- 
nian with three-body interaction ! 2 * 3 ! 4 ! These topological 



phases admit anyonic excitations with non-abelian braid- 
ing statistic. Of special interest are also spin systems and 
bosonic Hamiltonians with complex many-body interac- 
tions, such as ring exchange model , which are expected 
to give rise to exotic phasesP^E] Three-body interac- 
tions are also an essential ingredient for systems wi th a 
low energy degeneracy characterized by string netsj^^ 
which play an important role in models for non-abelian 
topological phases. The main challenge is then to identify 
experimental accessible systems, where the two-particle 
interaction V(r) can be controlled, independent of the 
three-body interaction W(ri, Tj ■, rfc). 
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FIG. 1: Strengths of the dominant three-body interactions 
Wijk appearing in the Hubbard model for different lattice 
geometries: (a) one-dimensional setup (b) two-dimensional 
square lattice (c) two-dimensional honey comb lattice. The 
characteristic energy scale Wo = 72_D.R0/ffl 6 is discussed fol- 
lowing Eq. 1 14 1. 



In the present work we analyze the effective interac- 
tion potential in a many-body system of polar molecules. 
Recently, systems of polar molecules in the rovibra- 
tional ground state have attracted a lot of interest re- 
cently due to their rich internal structure and the pres- 
ence of large permanent dipolc moments, which give rise 
to dipole-dipole interactions and offer the possibility to 
tune the inte raction with s tatic electric fields and mi- 
crowave fields EEMMiElESl The techniques for trapping 
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and cooling of polar molecules with the goal to cre- 
ate quantum degenerate ground state molecules are cur- 
rently developed in several laboratories! 17 1 18 1 19 1 20 1 21 1 22 1 ^ 
We show below that for an appropriate choice of static 
electric and microwave fields the effective interaction re- 
duces to the form as in Eq. (JT|) with tunable two-body 
and the three-body interactions. 

In particular, for polar molecules moving in an optical 
lattice we obtain the Hubbard model 

H = -jJ2b\b j +Y,-^-n i n j + -^ni»i»fc- (2) 

Here bi (b\) are destruction (creation) operators for a 
molecule on lattice site i, satisfying canonical commuta- 
tion (anti-commutation) relations for bosonic (fermionic) 
molecules, and a hard-core on-site repulsion is implied for 
bosons. The first term in Eq. ([2| is a hopping term (ki- 
netic energy) , while the last two terms describe two-body 
and three-body interaction terms (n, = b\bi). The dif- 
ferent strengths of the three-body interaction terms are 
illustrated in Fig. [T| We emphasize that our derivation 
of the Hubbard model Eq. ([2|, resulting in strong and 
tunable two and three-body interactions, will be based 
directly on the effective many-particle potential Eq. (JT|). 
This is in contrast to the common approach to derive ef- 
fective many-body terms from Hubbard models involving 
two-body interactions, which arc obtained in a J <C U 
perturbation theory, and are thus necessarily small. 24 
The main part of the present work is concerned with the 
microscopic derivation of the Hamiltonian in Eq. ([2| , and 
the tunability of the parameters by external fields. As 
an illustration, we analyze the simplest possible case of 
a one-dimensional Bose-Hubbard model with a dominant 
three-body interaction, and derive the phase diagram us- 
ing Bosonization techniques. 

I. EFFECTIVE INTERACTION POTENTIAL 

The internal structure of polar molecules with a closed 
shell electronic structure 1 E is given by the rotational 
degree of freedom, and its low energy excitations are well 
described by the rigid rotor Hamiltonian 

Hj& = B3i (3) 

with B the rotational constant, and the dimensionless 
angular momentum operator. Here, we are mainly in- 
terested in the ground state |0, 0),, and the first excited 
state manifold \l,m z )i (\j,m z )i denote the eigenstates 
of the rotor). The rotational levels couple to external 

electric fields via H£l — — d;E(i) with the dipole op- 
erator. Under a static electric field E = Ee z along the 
z-axis the degeneracy of the excited state manifold with 
j = 1 is lifted, see Fig. [2] In the following we denote 
the new dressed states by 1 0) ^ — > \g)i for the ground 
state, \l,0)i — > \e)i for the state with m z — 0, and 



|1, ±l)j — » \e±)i for the remaining two degenerate states, 
and the corresponding energies E g , E e , and E e ^±. 




FIG. 2: (a) Internal excitation energies for a single polar 
molecule in a static electric field E = Ee z . (b) Level struc- 
ture for Ed/B = 3: the circular polarized microwave field 
couples the ground state \g) with the excited state |e+) with 
Rabi frequency Q/h and detuning A. Applying a second mi- 
crowave field with opposite polarization (red arrow) allows us 
to lift the degeneracy in the first excited manifold (j = 1) by 
resonantly couple the state |e_) to the manifold with j — 2. 

We focus on a setup with a circular polarized mi- 
crowave field along the z-axis, which couples the states 
\g)i and |e+)$ with detuning A and Rabi frequency Cl/h. 
We assume that the degeneracy of the states |e_) and 
|e_|_) is lifted, e.g., by an additional microwave field cou- 
pling the state |e_) near-resonantly to the next state 
manifold inducing an AC shift, see Fig. [2] the degen- 
erate situation is presented in Appendix A. The internal 
structure of a single polar molecule is then described as 
a spin-1/2 particle via the identification of the state \g)i 
(|e+)j) as eigenstates of the spin operator Sf with posi- 
tive (negative) eigenvalue. In the rotating frame and ap- 
plying the rotating wave approximation, the Hamiltonian 
describing the internal dynamics of the polar molecule re- 
duces to 

(£-")-"* w 

with the effective magnetic field h = (f2, 0, A) and the 
spin operator Si = (Sf , Sf , Sf). The eigenstates of this 
Hamiltonian are denoted as = a\g)i + /3|e_|_), and 
= -f3\g)i + a\e+)i with energies ±VA 2 + f2 2 /2. In 
the following, we consider a blue detuned microwave field 
with A > 0, and molecule prepared into the state 
The preparation in this state is obtained by adiabatically 
turning on the microwave field for molecules initially in 
their ground state \g)i- 

The interaction between the polar molecules is deter- 
mined by the dipole-dipole interaction 

with Vij — Yi—Yj the separation between the particles. In 
the interesting regime with |ry| » (D/B) 1 ^ 3 , the dipole- 
dipole interaction restricted to the internal states |e+)j 
and \g)i expressed in the rotating frame, and applying the 
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rotating wave approximation reduces to H d — iJ° x + H* 
The first term takes the form (Sf 1 = Sf ± iSf) 



H ™ = ~lH yK r ;j) [s+sr + s+s-] 



(6) 



with the dipole coupling D — |(g|di|e+)| 2 , while the in- 
duced dipole moments d g — dEE g and d e — dEE &i+ gives 
rise to a second term 

H ? = \T, DU ( F « ) frgPi + VeQi] {VgP 3 + VeQj] ■ (7) 

Here, P, = 1/2 + Sf and Qi = 1/2 — Sf are the projectors 
on the ground and excited states, while r\ g — d g /yD and 
r/ e = d e j\[D are the dipole couplings. The anisotropic 
behavior of the dipole-dipole interaction is accounted for 
by v(v) = (1 — 3cos 2 #)/r 3 with 9 the angle between r 
and the z-axis. 

Next, we are interested in the effective interaction be- 
tween the polar molecules with each molecule prepared in 
the state Within the Born-Oppenheimer approxi- 

mation, we determine the eigenenergies of the internal 
Hamiltonian J2i Hq + Hd f° r fixed particle positions 
{r^}, and obtain the energy shift of the state adiabat- 
ically connected to the state \G) — of the non- 

interacting system. This energy shift is driven by the 
dipole-dipole interaction H d and strongly depends on the 
positions of the particles {r^} and therefore describes the 
effective interaction V t>e ({ri}). In the following, we de- 
rive this energy shift using perturbation theory in the 
dipole-dipole interaction H d , and derive the effective in- 
teraction potentials. The small parameter controlling the 
perturbative expansion is D/(a 3 |h|) = (R^/a) 3 with a 
the characteristic length scale of the interparticle separa- 
tion, |h| = \/ A 2 + ft 2 the strength of the effective mag- 
netic field, and the length scale R = (D/VA 2 + fi 2 ) 1 / 3 . 
In first order perturbation theory, we obtain the energy 
correction 

E W ({n})=\ [(a\+(3 2 Ve ) 2 ~a 2 p 2 }j2Dv(r tJ ), (8) 

which describes a dipole-dipole interaction between the 
particles. In addition, the energy shift in second order 
perturbation theory reduces to 

^({r,}) = J2 ^J^^ M^Mr^) 



y } N]2 [Dv{v l3 )\ 2 

^ 2%/a 2 + n 2 1 v J ' s 



K 3 



(9) 



and gives rise to a correction to the two-particle interac- 
tion potential and an additional three-body interaction. 
The matrix elements M and N take the form 

M = a(3 [(a 2 7y 3 + fi 1 ^) ( Ve - Vg ) + (/3 2 - a 2 )/2] , 



Therefore, the effective interaction potential V a!! up to 
second order in (R /a) 3 reduces to the form in Eq. |l]) 
with the two-particle interaction potential 



(10) 



N 



a z f3 2 (r) e -r) g y + l 



V(r) = ^Dv^ + ^DR 3 [u(r)} 2 , 

and the three-body interaction 

W(r u r 2 ,r 3 ) = l2 R 3 D [v{r X2 )v{r 13 ] (11) 

+Kr 12 )Kr 23 ) + 13)^23)] • 

The dimcnsionlcss coupling parameters are Xi = 
(a 2 Vg + /3 2 Ve ) 2 - a 2 /3 2 , A 2 = 2|M| 2 + |7V| 2 /2, and 72 = 
2|M| 2 . These parameters can be tuned via the strength 
of the electric field Ed/B and the ratio between the Rabi 
frequency and the detuning, fl/A, see Fig. [3] Of spe- 
cial interest are the values of the external fields, where 
Ai = 0, i.e., the leading two-particle interaction vanishes. 
Then, the interaction is dominated by the second order 
contribution with A2 and 72, which includes the three- 
body interaction, see Fig. [SJi, while small deviation away 
from the line Ai = allows us to change the character of 
the two-particle interaction. Note, that a n-body inter- 
action term (n > 4) appears in (n — l)-th order perturba- 
tion theory in the small parameter (Ro/a) 3 . Therefore, 
the contribution of these terms is suppressed and can be 
safely ignored in our context. 

The above analysis for \G) = Hi\+) (A > 0) provides 
a positive energy shift in second order £^ 2 )({ri}), as \G) 
corresponds to the highest energy state, and provides a 
repulsive interaction with A2 > and 72 > 0. In turn, in 
an analogous analysis for the lowest energy state IT|— ), 
is negative, which yields a change of sign of the 
coupling parameters A2 and 72. 

The validity of the effective interaction I/ ff ({rj}) is re- 
stricted to interparticle distances |r; — ry| > Rq. Here, 
we prevent particles to approach each other on shorter 
distances — Tj\ < i? by focusing on setups where the 
combination of interparticle potential for |r^ — r ^ | > Rq 
and the transverse trapping potential produces a strong 
repulsive barrier. For sufficiently strong barrier height, 
thermal activation across the barrier and quantum me- 
chanical tunneling through the barrier are then sup- 
pressed, and the particles are confined in parameter 
space to the region |r.j — rj\ > Rq. A setup provid- 
ing such a barrier is obtained by confining the parti- 
cles into two-dimensions by a strong transverse trapping 
potential along the z-axis with transverse trapping fre- 
quency lu± = h/raa 2 ^ as provided for example by an 
optical lattice. The condition for an efficient barrier 
in this two-dimensional setup has been recently worked 
out for po larize d molecules with leading dipole-dipole 
interaction! 25 * 26 ^ However, the analysis can be generalized 
to the present situation with the interaction potential 
Kft({ri}), if the two-particle potential V(r) is sufficient 
repulsive, i.e., Ai > — A2(i?o/a) 3 - For f&X > D/R 3 , 
an estimate of the rate for two-particles to penetrate 
the barrier up to a distance jr, — r,-| < R is pro- 
vided by r ~ (h/ma 2 ) exp(—2S E /h) with the semiclas- 
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FIG. 3: (a)-(c) Strength of the interaction parameters Ai, A2, 
and 72 as a function of the external fields Ed/B and Q/A. 
The leading dipole-dipole interaction vanishes for Ai = 
(dashed line in (b) and (c)), and the second order contri- 
butions dominate the interaction, (d) Strength of A2 (dashed 
line) and 72 (solid line) along the line in parameter space with 
Ai = 0. 



sical action S E /h ~ \f Dm / 'RqH 2 . This exponential sup- 
pression confines particle in parameter space to the re- 
gion |r j — Tj\ > Rq for realistic parameters with polar 
molecules, see below. 

The low-energy many-body theory now follows by com- 
bining the kinetic energy of the polar molecules with the 
effective interaction V ott within the Born-Oppcnhcimcr 
approximation and the external trapping potentials V T 



trapping potential. We describe the lattice structure 
with lattice spacing a by a set of vectors {Ri} account- 
ing for the minima of the optical lattice; depending on 
the optical lattice, we can generate one-dimensional and 
two-dimensional systems. The mapping to the Hubbard 
model follows the standard procedure!^ (i) Solving the 
single particle problem in the presence of the optical lat- 
tice provides the Wannier functions w(r) for the lowest 
Bloch band and determines the hopping energy J. The 
Wannier functions describe localized wave functions with 
characteristic size ao- (ii) We express the field operator 
ip in terms of the Wannier functions in the lowest band 
V-^(r) = J2i w ( r — RiWi with b\ the creation (annihi- 
lation) operator at the lattice site R^. In order to simplify 
the discussion, we consider bosonic particles satisfying 
bosonic commutation relations. Expressing the Hamil- 
tonian Eq. ( [T2| in second quantization and inserting the 
bosonic field operator ip, maps the system to the Bose- 
Hubbard model. However, the on-site interaction, which 
derives for cold atomic gases from the pseudopotential, 
requires a special discussion in the present situation. In 
the experimentally interesting regime, we have the fol- 
lowing separation of the length scales: ao < Ro < a. 
With the above discussion, that the parameter space of 
the system is confined to interparticle distances larger 
than i?o, this implies that if a particle is at site Rj, 
the hopping rate for a second particle to tunnel to this 
site is strongly suppressed. As the initial system has no 
doubly occupied sites, a convenient way to express this 
conditional hopping is to describe the bosons as hard- 
core bosons. Consequently, no on-site interaction term 
is present, and we obtain the Bose-Hubbard model in 
Eq. ([2]). The interaction parameters E/y and Vijk derive 
from the effective interaction V*({ri}), and in the limit 
of well- localized Wannier functions (do/a <C 1) reduce 
to U i:i = ViR, - Rj) and W ljk = W(Ri,Rj,R k ). The 
decay of these interactions with interparticle separation 
takes the form 



(12) 



Note, that the Hamiltonian is independent of the statis- 
tics of the particles and therefore, valid for bosonic and 
fcrmionic polar molecules. In the strongly interacting 
regime, where the interaction energy dominates over the 
kinetic energy, it is expected that the ground state of the 
many-body system is determined by crystalline phases.^ 



II. HUBBARD MODEL 

Applying an optical lattice provides a periodic struc- 
ture for the polar molecules described by the Hamilto- 
nian Eq. (12) and allows us to derive Hubbard models 



with unconventional and strong nearest-neighbor inter- 
actions. We focus on the above setup, where the sta- 
bility of the system is obtained by a strong transverse 
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W ijk = Wo 
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Rj| 3 |Rj 
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perm 



(13) 



(14) 



Here, U = AiL>/a 3 , U x = X 2 DRl/a 6 , and W = 
72-D-Rg/a 6 denote characteristic energy scales. The domi- 
nant contributions and strengths of the three-body terms 
in different lattice geometries are shown in Fig. [l] 

In the following, we estimate these energy scales for 
LiCs with a permanent dipole moment d = 6.3Debye. 
Assuming an optical lattice with lattice spacing a ~ 
500nm, the characteristic dimensionless parameter deter- 
mining the ratio between the interaction energy (E int = 
D/a 3 ) and the characteristic kinetic energy within the 
lattice (E kin = h 2 /ma 2 proportional to the recoil energy) 
becomes — Dm/h 2 a ss 55. The leading dipole-dipole 
interaction can gives rise to very strong nearest-neighbor 
interactions with Uq ~ 55£" kin . On the other hand, tun- 
ing the parameters via the external fields to Ai =0 the 
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TABLE I: Parameters of the sine-Gordon term 



characteristic energy scale for the three-body interaction 
becomes W « (R /a) 3 E kin . Then, controlling the hop- 
ping energy J via the strength of the optical lattice allows 
us to enter the regime with dominant the three-body in- 
teractions. 



III. THREE-BODY INTERACTIONS IN ID 

As an illustration, we study the one-dimensional sys- 
tem with leading three-body interaction. Tuning the ex- 
ternal fields to a point with U\ = —Uq, the Hamiltonian 
takes the form 



H 



-J 



+i 



fct 



■W)^ni-inin i+ i, (15) 



where we have kept only the nearest-neighbor terms. The 
additional interaction terms play a minor role in the fol- 
lowing analysis. Using a Jordan- Wigner transformation 
allows us to map the hard-core bosons onto an interact- 
ing Fermi syste m, wh ich can be studied using Bosoniza- 
tion techniques) 28 * 29 ! the details of the calculation are 
presented in Appendix B. As a final result, we obtain 
that the low energy Hamiltonian is described by the sine- 
Gordon model for the particle densities n = 1/3, n = 1/2, 
and n — 2/3, 



H 



hv 
Y 



dx KTL 2 + — (d x <S>) 2 +u / dx^-^ cos (/?$) . 
K J J TT z a z 

(16) 

The parameters u and (3 for the sine term are given 
in Table [I] The renormalized Fermi velocity is v = 
2Jasin(j:n)/hK, while the dimensionless Luttinger pa- 
rameter takes the form K = 1/ \f\ + k with 



W 2n 3 — 2 cos(27rn) — cos(47rn) 



2J 7T 



sin(7m) 



(17) 



The Hamiltonian Eq. (16) allows now for the derivation 
of the expected phases diagram. 

For weak interactions W/2J < 1, the three-body 
interaction shifts the hard-core bosons away from the 
Tonks gas limit (K = 1) into the correlated regime 
with K < 1. This regime is characterized by alge- 
braic correlation functions, e.g., the off-diagonal correla- 
tion between the bosons decays as (6j&t) ~ \i — j\~ 1 ^ 2K , 
and the density-density correlation ((n, — n)(nj — n)) ~ 

cos[27m(i — — j\ 2K . 30 At three values of the den- 

sity n — 2/3, n — 1/2, and n — 1/3, the sine-Gordon 
term drives an instability towards a gapped phase, see 



Fig. [4] The critical value K for this instability is given 
by K = 8tt//3 2 . In the following, we use this criterion 
to identify the instabilities towards three solid phases for 
W/2J ~ 1; the exact determination of the critical inter- 
action strength requires an analysis beyond the scope of 
this paper. These solid phases are characterized by a bro- 
ken translational symmetry, where the following observ- 
ables allow us to distinguish the different ground states 



cos 27r jn + V 4:t:^i , (18) 
(2nj + ir)n + V4:Tr$i (19) 



cos 



with <1>; characterizing the different ground states. The 
first observable describes a density wave with wave vec- 
tor k = 27rn/a present in conventional solids, while the 
second observable accounts for a bond order with wave 
vector k — 2im/a, see Fig. [4Jd for an illustration of the 
different order appearing in the ground state. 




density wave k — 2ir/3a 



bond order k = ir/a 



bond order k — 2n/3a 



FIG. 4: (a) Sketch of the W/ J-n phase diagram: The dashed 
contours correspond to constant values of K. The first insta- 
bility at n — 1/2 appears for K — 1/2 (thin solid line), while 
the instabilities at n = 1/3 and n = 2/3 appear at lower 
values of K. Note that the interaction strength approaches 
W/2J ta 1 at the position of the instabilities, and we expect 
the exact position of the instability to be strongly renormal- 
ized. (b) Illustration of the three different solid order char- 
acterizing the gapped phases: density order with wave vector 
k = 2n/3a, bond order with k = n/a, and bond order with 
wave vector k — 2ir/3a. 

At half filling n— 1/2, the instability towards a gapped 
phase appears for K = 1/2. Within the gapped phase, 
the sine-Gordon term in Eq. ( 16 1 determines the long- 



distance behavior, and the phase field <f> is predominantly 
pinned within a minimum of the cos(\/l67r<i>) term. The 
minima take the form $>i = ir(2l + l)/\/167r, and char- 
acterize the different ground states in the gapped phase. 
These ground state can be distinguished by the bond ob- 
servable in Eq. (19 1. Therefore, we obtain a two- fold 



degenerate phase with a broken translational symmetry: 
the bond correlation function exhibits a long range or- 
der at the wave vector k — n/a, while the density (nj) 
remains uniform in this phase. The corresponding phase 
in spin systems is denoted as a spin-Peierls phaseP^l 
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In turn, for the densities n = 1/3 and n = 2/3 the in- 
stability appears at lower values of K . The ground states 
of the gapped phase are then characterized by the minima 
of cos(v / 3&7r < i ) ), which are given by = tt{21 + l)/y/36n. 
The different ground states can be distinguished by the 
density and bond observable providing a three-fold de- 
generate phase with a density wave and a bond order at 
the wave vector k — 2it /3. The appearance of a density 
wave with k — 2tt/3 is a special property of the three- 
body interaction. It can be well understood for n = 2/3 
in the limit W/2J ^S> 1: then the ground state takes the 
form Ilj&g^g^jjO) and describes a perfect density wave. 

Within the above Bosonization approach, the addi- 
tional interaction terms beyond nearest-neighbor three- 
body interaction only provide a small renormalization of 
the coupling parameters in Eq. (16 1, and therefore play 



a minor role in the qualitative discussion of the insta- 
bilities. However, in the limit W/2J ^> 1 they become 
important and can give rise to additional solid phases. 



IV. CONCLUSIONS AND OUTLOOK 

A many-body system of cold polar molecules, whose 
rotational states are dressed by an external static electric 
field and microwave field, can be described by an effec- 
tive Hamiltonian, where three-body interactions play the 
dominant role. The derivation of this effective Hamil- 
tonian involving interaction potentials as a series of 
two-body, three-body etc. interactions is based on a 
microscopic understanding of the underlying molecular 
physics, and is systematic in the sense that the series ex- 
pansion is well controlled. A unique feature of the system 
is that, as a function of the external fields, the two-body 
interactions can be tuned from repulsive to attractive, 
and even switched off, while the three-body terms remain 
repulsive and strong. For molecules trapped in an optical 
lattice this leads to Hubbard models with tunable near- 
est neighbor two-body interactions and repulsive three- 
body terms. Models of this type have appeared in the 
recent discussion of exotic quantum phases, in particular 
in the context of topological quantum phases and quan- 
tum computing, and we see molecular quantum gases as 
a realistic experimental route, which provides the basic 
building blocks and techniques towards the study of these 
phenomena. 



APPENDIX A: DEGENERATE STATES 

For a setup with |e + ); and |e_)i degenerate, it is neces- 
sary to keep the three states and |e_)j for the 
perturbative calculation of the Born-Oppcnhcimcr po- 
tentials. The leading contribution £i({r.;}) is not modi- 
fied, while the following term in the dipole-dipole inter- 
action provides a non-vanishing contribution in second 



order perturbation theory, 
A#r = \ £ D M T *N T ? S 7 + M(r 4 ,)?TS+] (Al) 



with the operators = |e_)(g| and = |<?)(e_| cou- 
pling the ground state with the excited state and the 
potential /Lt(r) = —3 (x — iy) 2 /2r 5 . We therefore obtain 
a correction to the two-body interaction potential 



(A2) 



AV(r) = X 3 DR 3 | M (r)| 2 , 
and the three-body interaction 

n e>3 3 

W(r 1 ,r 2 ,r 3 )=73-P £ [mM Vfo*)- (A3) 
The dimcnsionless parameters A3 and 73 take the form 



ofl4 2 V^ 2 + A 2 4fl2 V^ 2 + A 2 

A 3 = 2f3 4 a 2 - — — — + 4a 4 /? 2 



73 = 4a 4 /3 



A + 3Vtt 2 + A 2 
,., + A 2 



A + VA 2 + ft 2 



A + x/A 2 + ft 2 ' 
and depend on the external fields Ed/B and il/A. 

APPENDIX B: BOSONIZATION 

Using the equivalence between hard-core bosons and 
a spin-1/2 systems, allows us to map the Hamiltonian 
in Eq. ( 15 1 to a fermionic model with Fermi operators 
Cj (cj) via a Jordan- Wigner transformation. Follow- 
ing the standard Bosonization procedure, we express the 

fermioni c oper ators via slowly varying left and right mov- 
ing fields™] 



d ~ Va [e~ ikFXi R(x) + e lk ^L{x)] 



(Bl) 



with k F = nn I a the Fermi momentum and n the averaged 
particle density. Here, Xi describes the position of the 
lattice site i, while a denotes the lattice spacing. The 
fields R(x) and L(x) are slowly varying and smooth on 
distances a in the continuous variable x ~ Xj with 



(B2) 
(B3) 



R(x) — exp ( i V 47r</> 

'2ira 



L(x) 




exp i^/47T(j)j 



Then, the Hamiltonian for non-interacting fermions maps 
to the Luttinger liquid Hamiltonian 



IT 



J dx[ll 2 + (d x <S>) 2 ] (B4) 



with the bosonic field <1> = 4> + 4> and n the momentum 
conjugate satisfying the canonical commutation relation 



7 



[<&(x),n(x')] = iS(x - x'). Here, v F = 2Jsm(irn)/h de- 
notes the Fermi velocity of the non-interacting Fermi sys- 
tem. In order to study the influence of the three-body 
interaction, we split the density operator into its mean 
value and the fluctuations, i.e., rii = n + Arii. Then, the 
interaction Hamiltonian reduces to H int = H± + H2 + H 3 
with 

H 2 = Wn,y) [2An»Ant + i + Ani_iAn i+ i] , (B5) 

i 

H 3 = Wj2 An i-l An i An i+l> ( B6 ) 

i 

while Hi denotes a shift in the chemical potential and 
can be dropped. The density operator expressed in the 
right and left moving fields takes the form 

An, = a [d x $ + e 2ik * Xi M + (x) + e~ 2ik * Xi M~ (x)] 

with 

M+(x) =: R\x)L{x) :, M~(x) =: h\x)R{x) : . 

Expressing the interaction Hamiltonian H int in terms of 
the bosonic fields <&, special care has to be taken due to 
the normal ordering : : of the operators. Furthermore, 
we note that some terms exhibit a fast oscillations, which 
can be ignored in the low energy Hamiltonian. Then, the 
interaction Hamiltonian H2 takes the form 

H int = Wna I dx- [3 - 2cos(2?m) - cos (47m)] (d x $) 2 . 



Furthermore, at half filling with n = 1/2 an additional 
term appears due to cancellation of different oscillating 
terms, which takes the sine-Gordon form 



H 



sG 



dx^—r cos(/3$) 



(B7) 



with u = Wna and [3 — v 16-7T. The interaction H 3 only 
contributes less relevant terms, except at fillings n = 1/3 
and n = 2/3. Then an additional terms appears, which 



takes the sine-Gordon form in Eq. (B7l with parameters 
u = Wa/ir and /3 = V367T. 
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